Density functional theory and demixing of binary hard rod-polymer mixtures 
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A density functional theory for a mixture of hard rods and polymers modeled as chains built of 
hard tangent spheres is proposed by combining the functional due to Yu and Wu for the polymer 
mixtures [J. Chem. Phys. 117, 2368 (2002)] with the Schmidt's functional [Phys. Rev. E 63, 50201 
(2001)] for rod-sphere mixtures. As a simple application of the functional, the demixing transition 
into polymer-rich and rod-rich phases is examined. When the chain length increases, the phase 
boundary broadens and the critical packing fraction decreases. The shift of the critical point of a 
demixing transition is most noticeable for short chains. 



The addition of the nonadsorbing polymer to a mono- 
disperse suspension of colloidal particles can lead to a 
phase separation due to depletion interactions aris- 
ing from a tendency of the system to reduce the vol- 
ume excluded to the polymer coils. One of the simplest 
theoretical models taking into account this phenomenon 
is the so-called Asakura-Oosawa (AO) model of colloid- 
polymer mixtures 0] where the ideal polymer coils (mod- 
eled as spheres) can freely interpenetrate each other but 
the polymer-colloid and colloid-colloid interactions are 
of the hard sphere type. Initial studies of such systems 
focused on the bulk phase behavior 0, , however re- 
cently developed density functional theory (DFT) for the 
AO model |5| initiated investigations of inhomogeneous 
colloid-polymer mixtures. When brought close to a hard 
wall, such mixtures may develop a sequence of layering 
transitions in the partial wetting regime prior to a tran- 
sition to complete wetting [(J Q • 

Similar mechanism of fluid-fluid phase separation can 
be found if other mesoscopic particles such as hard rods 
are used as depletant agents. Bolhuis and Frenkel (BF) 
|8| used computer simulations and free volume theory 
[4| to study bulk phase behavior of mixtures of colloidal 
hard spheres and vanishingly thin hard rods. They found 
a surprisingly good agreement (c.f. Fig. 3 in Q) between 
the free volume theory and Gibbs ensemble Monte Carlo 
simulations. Although the vanishing thickness and vol- 
ume constitutes a significant simplification (e.g. it rules 
out the isotropic-nematic transition), the BF model is 
relevant for some experimental systems e.g. mixtures of 
silica spheres and silica coated bohemite rods 0, . In 
order to describe inhomogeneous hard rod-hard sphere 
systems Schmidt ^lj proposed a density functional the- 
ory which in its structure closely resembles the AO func- 
tional .5]. The functional for the BF model incorpo- 
rates the exact low-density limit and yields the same 
equation of state as in Ref. Moreover, entropic sur- 
face phase transitions found previously in model colloid- 
polymer mixtures close to a hard wall were also recently 
encountered in hard rod- hard sphere mixtures 01 ■ This 



further demonstrates deep similarities between the two 
models. 

The aim of the present work is to construct a functional 
for a mixture of vanishingly thin hard rods and polymers. 
Such systems can be regarded as simple microscopic mod- 
els of the liquid crystal-polymer mixtures. The functional 
is constructed by combining the Schmidt's functional for 
the BF model with the Yu and Wu ^| functional for 
mixtures of polymeric fluids. To make this conjecture 
we take the advantage of the fact that both function- 
als underlie the fundamental measure theory (FMT) of 
Rosenfeld |T^|. As a simple application we investigate 
bulk phase diagrams resulting from the proposed theory. 

Consider a mixture of hard, vanishingly thin needles 
(species N) of length L and polymers (species P) mod- 
eled as chains composed from M tangentially bonded 
hard-sphere segments of diameter a. The hard-sphere 
monomers building up the chains are freely jointed i.e. 
they can adopt any configuration as long as it is free 
of the intermolecular and intramolecular overlap. The 
interaction potential between needles Vnn — for all 
separations, while the pair potential between a polymer 
segment and a hard rod, Vpn, and between two polymer 
segments, Vpp, is of a hard-core type i.e. is infinite if 
a pair of objects overlap and zero otherwise. The grand 
potential of such system as a functional of local densities 
of polymers pp(R) and needles pf^(r,u:) can be written 
as 

n[p P (R),p N (r,u>)} = Fid[pp(R),PN{r,oj)]+ 

F ex [p P (K),p N (r,u:)} + J dRpp (R) (R) - p P ) 



du) 



<h I -^PN(r,u){V e { xt\r,u) - p N ) , 



(1) 
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where R = (ri, r 2 , • • ■ , tm) denotes a set of coordinates 
describing the segment positions, uj describes the orienta- 
tion of the rod, V^'(R), pp, V^P (r, w) and pn are the 
external and the chemical potentials for polymers and 
rods, respectively. The ideal part of the free energy is 
known exactly 

l3E l d[pp(R),p N {r,u J )]= /3 /dRp F (R)Vb(R)+ 



2 



/dRp P (R)KApPp(R))-l] 



4~ 



pjv(r,a;)KA^iv(r,w))-l], (2) 



where Ap and A^v are the thermal wavelengths for poly- 
mers and needles, respectively. The total bonding po- 
tential Vb (R) is a sum of bonding potentials Vb between 

the segments Vf,(R) = Yli=i 1 ^b(l r i+i — r i|) and for 
the tangential hard spheres this contribution satisfies 
exp[-/3Vb(R)] = nf'7 1 5(|r. i+1 - r,| - ct)/4™ 2 . 

Within the framework of the fundamental measure 
theory the excess free energy density $ is expressed as 
a simple function of the weighted densities . We 
assume that <£> can be represented as a sum of the 
orientation-independent polymer contribution, &pol, 
and the orientation-dependent polymer-needle contribu- 
tion, $pat. For the orientation-independent contribution 
we use an extension of the Rosenfeld's FMT to the poly- 
meric fluids [l3i |. where the polymer excess free energy 
density $pol is assumed to be a functional of only seg- 
ment densities pps(r) defined as 



m m . 

p PS {v) = PPsAr) = E J dR ^ r ~ ^)PP( R ) 



(3) 



where pps.i(r) is the local density of the segment i of 
the polymer. Following [l3| we assume that $poi can 
be split into the hard-sphere &hs contribution resulting 
from the hard-sphere repulsion of polymer segments and 
$p contribution arising due to the chain connectivity. 
Taking into account the assumptions mentioned above 
we have $ = $paj + &hs + $p. The total excess free 
energy of the inhomogeneous system is obtained via the 
integration of the excess free energy density 

PF ex = JdvJ ^$ PiV ({n«})+$H5({n£ p )})+$p({n^)}). 

For the hard-sphere part we choose the original 
Rosenfeld expression 



e p; 
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*Hs({n ( a P) }) = -n^'\n(l-nD+ (5) 
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The polymer weighted densities n a (r) are evaluated via 
spatial convolutions 



i p \v)= dv' PPS {r')wi P \r-v') 



(6) 



( P) 

where the weight functions Wa (r) are given by 



4 F) (r) = 0(|-|r|) , wP(r) = 5(^-\r\),(7) 
«^(r) = - |r|) , -r ) W = ^ ! (8) 



W " W Ira* 



The contribution due to the chain connectivity is evalu- 
ated using the Wertheim's first-order perturbation the- 
ory for a bulk fluid 0] and extended (using FMT-style 
weighted densities) by Yu and Wu to inhomogeneous sys- 
tems 



<&p(R p) }) = 



1 - M 

M 



n P CHy hs {a,{n {P) })], (10) 



where Q = 1 — riy^ • n |/2 /^P) 2 1 while 2//i S is connected 
with the Carnahan-Starling expression for the contact 
value of the radial distribution function of a hard-sphere 
mixture 



y hs (aAn iP) }) 
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l-n {P > 4(1 -7i^ ; ) 2 72(1 -nf) 



(P) 



(11) 

We note here that both <£>#s and <£> p are independent on 
the local density of rods. 

To specify the polymer-needle contribution we make 
use of Schmidt's DFT for hard rod-sphere mixtures [ll| 
and write the excess free energy density due to vanish- 
ingly thin needles as 



(AT) (PAT) 

$PA r ({n«}) = -4 iV) Ml-4 P) )+ 1 " 2 



1 -71 



(P) 



(12) 



The equation above implies that the polymer-needle con- 
tribution to the excess free energy stems solely from 
the hard-core repulsion between the needle and the hard 
sphere forming the polymer segment. 

The needle weighted densities, Tia ; are obtained 
through spatial convolutions of the needle local den- 
sity and the corresponding orientation-dependent weight 
functions 



(13) 

while the "mixed" polymer segment-needle weighted den- 
sity, ri2 PN \ is obtained via spatial convolution of the 
polymer segment density and an orientation-dependent 
weight function 

4 PN \v,u;) = Jdv' p PS (v')wi PN \v-r',u;). (14) 

Corresponding weight functions are given as ^l| 

w[ N) (r,u) = \ ( L/2 dW(r + wl) , 

4 J-L/2 



w 1 N) (t,lj) = -[5(v + u>l) + 5(v-u>l)} , 



w {PN) (r,u>) = 2|u,£ ) (r)-"| 



(15) 



This completes the prescription for the functional. The 
present theory reduces to the Schmidt's functional ^lj if 
M = 1 and to Yu and Wu functional [l3| if the density 
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FIG. 1: Phase diagrams for hard rod-polymer mixtures. Solid 
and dashed lines denote the binodals and spinodals, respec- 
tively. Black circles indicate the critical points of a demixing 
transition. The main plot shows phase diagram for a mixture 
of hard rods with size ratio q = L/a = 1 and hard-sphere 
10-mers, in terms of polymer packing fraction np = ^a 3 Mpp 
and the scaled needle density p^L 2 a. Open diamonds denote 
the critical points for a mixture of hard-sphere 10-mers and 
hard rods with size ratios q = 20, 15, 8, 5 and 0.1 (from top to 
bottom). The inset shows the polymer packing fraction-scaled 
needle reservoir density representation of the phase diagrams 
for a mixture of hard rods with size ratio q = 1 and polymers 
with M = 1, 10 and 100 (from top to bottom). 



of rods pn — 0. Note also that, similarly to |Tl|, the 
functional is linear in the local density of rods. 

Although the proposed DFT is intrinsically designed to 
study inhomogeneous systems, here we restrict ourselves 
to the isotropic bulk phases where the density of both 
species are constant. In this case the scalar weighted den- 
sities become proportional to the bulk densities of needles 
and polymers, while the vector weighted densities vanish. 

Consequently from Eq. we obtain nj, 
d/Vs, with^ =*/6o* 

1 . Likewise Eq. ^5] yields 

= 1. Finally Eq. ITU leads to 
^(pjv) _ ^(p) p ps _ After inserting the above expressions 
for the weighted densities into the Eq. and Eqs. 1101121 
the total free energy per unit volume, is evaluated 
as $ v = <f> H s + $p + ®pn + f3~ 1 pp[ln(p P A 3 p ) - 1] + 
PJv[ m (/°A'A 3 v ) — !]■ From <&„ the pressure, p, and the 
chemical potentials of both species are easily calculated 



and £o 

where £^ = L/A and £q N 



& P) Mpp = 
= a/2 
Pn 



Pp = 



j=P,N 1 3 



0Hj = 



(16) 



hard rods undergoes entropically driven demixing transi- 
tion to polymer-rich (rod-poor) and polymer-poor (rod- 
rich) phases. The coexsisting equilibrium densities (bin- 
odals) were obtained by solving simultaneously equa- 
tions for the equality of pressures and chemical poten- 
tials in two phases. The spinodals delimiting the regions 
stable against fluctuations of density and composition 
were evaluated from the condition det[<9 2 <fr v /dpidpj] = 0, 
i,j = P, N. The critical points were evaluated from 



dpp 
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N 



N 
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d 2 $ v 
dp% 



0. 



(17) 



Under appropriate conditions a mixture of polymers and 



The above equation arises from the fact that at the crit- 
ical point the tie-line connecting the coexisting densities 
becomes tangent to the spinodal line [r7|. 

In Fig. n we show examples of binodals (solid lines), 
spinodals (dashed lines) and critical points (black circles) 
resulting from the present theory. The phase diagrams 
plotted in the polymer packing fraction, rjp — ^a 3 Mpp 

versus scaled needle reservoir density, p^ L 2 a represen- 
tation (see the inset) were evaluated for systems with 
M = 1, 10, 100 for constant q = L/a = 1. For the special 
case M = 1 (the uppermost diagram) the system reduces 
to the BF model jlq . 

As the chain length is increased, the phase boundary 
becomes more asymmetric and the critical point moves 
towards lower polymer packing fractions and towards 
lower reservoir needle densities. Another interesting fea- 
ture observed in |ll| is that when the polymer pack- 
ing fraction- actual needle density representation instead 
of the polymer packing fraction- reservoir needle density 
representation is chosen, the critical points for different 
size ratios q are located on the spinodal of the system 
with q = 1. We find the same behavior also in the present 
model i.e. for M > 1. In the main plot of Fig. ^ we show 
phase diagram for a mixture of 10-mers and hard rods. 
The critical points for mixtures with different size ratios 
(open diamonds, from top to bottom for q = 20, 15, 8, 
5 and 0.1) lie on the spinodal for the system with q = 1 
(dashed line). 

Let us now consider the limit of very long polymer 
chains. Intuitively one could argue, that in this regime 
the phase behavior of a needle-polymer mixture should 
barely depend on the rod length to polymer-segment di- 
ameter ratio q = L/a because the physical dimensions of 
the polymer e.g. the gyration radius, R g become much 
bigger than the rod elongation. This scenario is captured 
within the present theory. In Fig. [2] we show the phase 
diagram for a mixture of hard-sphere monomers (M = 1) 
and hard rods for the size ratio q = 1. The critical points 
for different size ratios q (black diamonds) are located on 
the spinodal. Open diamonds denote the critical points 
for the systems with M = 2, 5, 20, 100, and 10000, 
respectively. As the chain length increases the critical 
points for the systems with different size ratios but the 
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FIG. 2: The phase diagram for a hard rod-hard sphere 
monomer (M — 1) mixture with the size ratio q = 1. Solid 
and dashed line denote the binodal and the spinodal, respec- 
tively. Black diamonds on the spinodal indicate the critical 
points of a demixing transition for monomer-hard rod mix- 
tures with different size ratios q = 8, 5, 1 and 0.4 (going from 
top to bottom). Open diamonds denote the critical points of 
a demixing transition for a mixture of polymers consisting of 
M monomers and hard rods for fixed q and for M = 2, 5, 
20, 100 and 10000 (starting from the right hand side). Dot- 
ted lines serve as a guide to the eye and connect the critical 
points for the systems with the same size ratio q. 



same M come closer together and for M — 10000 they 
virtually merge (on the figure scale) into one point. In 
the limit M — > oo the critical scaled needle density tends 
to the value 2.5526 for all q while the critical polymer 
packing fraction tends to zero. Similar limiting behavior 
was found in mixtures of spherical colloids and polymers 
with excluded volume interactions 0, 0, ■ 

In conclusion, in this work we propose a density func- 
tional theory for a mixture of vanishingly thin hard rods 
and polymers modeled as chains built of hard tangent 
spheres. The functional is constructed by combining the 
functional due to Yu and Wu for polymer mixtures fl3| 
with Schmidt's functional for hard rod-sphere mixtures 
|ll| . The proposed theory predicts a demixing transition 
similar in its nature to that observed for sphere-rod sys- 
tems. The present functional is well suited to study inho- 
mogeneous systems. It would be of interest to consider a 
fluid-fluid interface or to investigate surface phase transi- 
tions such as entropic wetting or layering that have been 
discovered in colloid-polymer and colloid-rod mixtures. 
It is also straightforward to incorporate the Onsager limit 
[2ll| of the needle contribution to the functional generat- 
ing thus a more sophisticated model of inhomogeneous 
liquid crystal-polymer mixtures. Some of these topics 
are already under study in our laboratory. 
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